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oo , Abstract 
(N 

I The fully consistent relativistic proton-neutron quasiparticle random phase approximation (PN- 

RQRPA) is employed in the calculation of /3-decay half-lives of neutron-rich nuclei in the N~50 
and N~82 regions. A new density-dependent effective interaction, with an enhanced value of the 
nucleon effective mass, is used in relativistic Hartree-Bogoliubov calculation of nuclear ground 
^ I states and in the particle-hole channel of the PN-RQRPA. The finite range Gogny DIS interaction 



o 



is employed in the T = 1 pairing channel, and the model also includes a proton-neutron particle- 



^ . particle interaction. The theoretical half-lives reproduce the experimental data for the Fe, Zn, Cd, 

\ and Te isotopic chains, but overestimate the lifetimes of Ni isotopes and predict a stable ^^^Sn. 

PACS numbers: 21.30.Fe, 21.60.Jz, 23.40.Hc, 26.30.+k 
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I. INTRODUCTION 



The latest theoretical and computational advances in nuclear structure modeling have also 
had a strong impact on nuclear astrophysics. More and more often calculations of stellar 
nucleosynthesis, nuclear aspects of supernova collapse and explosion, and neutrino-induced 
reactions, are based on microscopic global predictions for the nuclear ingredients, rather 
than on oversimplified phenomenological approaches. The nuclear input for astrophysics 
calculations necessitates the properties of thousands of nuclei far from stability, including 
the characteristics of strong, electromagnetic and weak interaction processes. Most of these 
nuclei, especially on the neutron-rich side, are not accessible in experiments and, therefore, 
many nuclear astrophysics calculations crucially depend on accurate theoretical predictions 
for the nuclear masses, bulk properties, nuclear excitations, {n, 7) and (7, n) rates, a- and 
/3-decay half-lives, fission probabilities, electron and neutrino capture rates, etc. 

The path of the r-process nucleosynthesis runs through regions of very neutron-rich nuclei. 
/9-decays are particularly important because they generate elements with higher Z-values, 
and set the time scale of the r-process. Except for a few key nuclei, however, /3-decays of 
r-process nuclei have to be determined from nuclear models. One of the crucial questions 
for structure models is, therefore, a consistent microscopic calculation of /?-decay far from 
stability. Two microscopic approaches have been successfully applied in large-scale modeling 
of weak interaction rates of neutron-rich nuclei: the shell model and the proton-neutron 
quasiparticle random phase approximation (PN-QRPA). Shell model applications in nuclear 
astrophysics have recently been reviewed in Ref. jlj. In the calculation of /?-half lives, in 
particular, the principal advantage of the shell model is the ability to take into account 
the detailed structure of the /3-strength functions. In addition to large-scale shell model 
predictions for the half-lives of waiting-point nuclei at = 50, 80, 126 l|, the no-core shell 



model [2] and the shell model embedded in the continuum 3] have recently been applied to 
/9-decay of light nuclei. Large configuration spaces, however, prevent systematic applications 
of the shell model to heavy nuclei along the r-process path. 

When compared to the shell model, important advantages of an QRPA approach based 
on the microscopic self-consistent mean-field framework, include the use of global effective 
nuclear interactions and the treatment of arbitrarily heavy systems. Fully self-consistent 
QRPA calculations of weak interaction rates of neutron-rich nuclei, however, have only 
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recently been reported. In Ref. ^, in particular, a fully self-consistent QRPA has been 
formulated in Hartree-Fock-Bogoliubov (HFB) canonical single-particle basis. The canonical 
basis diagonalizes the density matrix and describes both the bound states and the positive- 
energy single-particle continuum. The formulation of the RQRPA in the canonical basis 
is particularly convenient because, in order to describe transitions to low-lying states in 
weakly bound nuclei, the two-quasiparticle configuration space must include states with 
both nucleons in the discrete bound levels, states with one nucleon in a bound level and 
one nucleon in the continuum, and also states with both nucleons in the continuum. The 
HFB+QRPA model based on Skyrme interactions, and with the inclusion of a finite-range 
residual proton-neutron particle-particle interaction, has been applied in the calculation of (3- 
decay rates for spherical neutron- rich r-process waiting-point nuclei ^. Based on a different 
energy-density functional, the self-consistent HFB plus continuum QRPA framework has 
been employed in a systematic calculation of the allowed and first-forbidden /3-decay rates 
for the r-process nuclei near = 50,80, 126 j^. It has been shown that the effect of the 
high-energy first-forbidden transitions is crucial for Z > 50, A^ ~ 82, and especially in the 
A^ = 126 region. 

In this work we report the first application of the relativistic PN-QRPA in the calcula- 
tion of /5~-decay rates for neutron-rich nuclei. Self-consistent relativistic mean-field models 
have been very successfully applied to nuclear structure, not only in nuclei along the valley 
of /3-stability, but also in exotic nuclei with extreme isospin values and close to the parti- 
cle drip lines. The relativistic mean-field framework has recently been extended to include 
medium dependent meson-nucleon vertices. One way to determine the functional form of 
the meson-nucleon vertices is from in-medium Dirac-Brueckner interactions, obtained from 
realistic free-space NN interactions. This microscopic method represents an ab initio de- 
scription of nuclear matter and finite nuclei. Another approach is purely phenomenologi- 
cal, with the density dependence for the a, u and p meson-nucleon couplings adjusted to 
properties of nuclear matter and a set of spherical nuclei. The resulting density-dependent 
meson-nucleon couplings are in qualitative agreement with those obtained in the microscopic 
approach. In Ref. Q we have extended the relativistic Hartree-Bogoliubov (RHB) model 
to include medium-dependent vertex functions. A phenomenological effective interaction, 
denoted DD-MEl, was adjusted to properties of nuclear matter and finite nuclei, and tested 
in the analysis of the equations of state for symmetric and asymmetric nuclear matter, of 
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ground-state properties of the Sn and Pb isotopic chains Q], and deformed nuclei It 
has been shown that, in comparison with standard non-hnear meson self-interactions, rel- 
ativistic Lagrangians with an explicit density dependence of the meson-nucleon couplings 
provide an improved description of asymmetric nuclear matter, neutron matter and nuclei 
far from stability. The relativistic random-phase approximation (RRPA), based on effective 
Lagrangians characterized by density- dependent meson-nucleon vertex functions, has been 
derived in Ref. In Ref. ^ the relativistic QRPA has been formulated in the canonical 
single-nucleon basis of the RHB model, and the corresponding proton-neutron relativistic 
QRPA (PN-RQRPA) has recently been employed in studies of isobaric analog resonances 
and Gamow- Teller resonances in spherical nuclei id. 

Our initial attempt to calculate (3~ decay rates of r-process nuclei by simply employing 
the DD-MEl interaction in the RHB plus PN-RQRPA framework was not successful. In 
general, the resulting half-lives were more than an order of magnitude longer than the 
empirical values. The reason is that DD-MEl, as well as most other successful relativistic 
mean-field interactions, has a relatively low effective nucleon mass, especially when compared 
to effective masses of Skyrme forces. This means that, because the density of states around 
the Fermi surface is low, in a self-consistent relativistic QRPA calculation of /5-decay the 
transition energies are low, and this result in long lifetimes. In some cases one could improve 
the results by adjusting the strength of the T = pairing interaction, which is not really 
known, but this is not the solution of the problem, and in any case it cannot be applied to 
doubly closed-shell nuclei. The problem of adjusting the effective nucleon mass in relativistic 
mean-field models is somewhat complicated, because this quantity is also related to the 
energy spacings between spin-orbit partner states. Since the latter are empirically well 
determined, the effective nucleon mass cannot be increased without extending the standard 
relativistic mean-field framework. In order to be able to calculate /3-decay half-lives, in this 
work we adjust a new relativistic mean- field density-dependent interaction with a higher 
value of the effective nucleon mass, by including an isoscalar tensor-coupling term in the 
model Lagrangian. In Sec. |n] we outline the theoretical framework and derive the new 
effective interaction DD-MEl*. In Sec. IIIII the relativistic PN-QRPA is employed in the 
calculation of /5-decay half-lives in the regions N^50 and N:^82. Sec. IIVI contains the 
summary, concluding remarks and an outlook of future applications. 
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II. OUTLINE OF THE THEORETICAL FRAMEWORK 



A. Proton-neutron relativistic quasiparticle random-phase approximation 

The relativistic QRPA has been derived in Ref . , in the canonical single- nucleon basis of 
the relativistic Hartree-Bogoliubov (RHB) model. The RHB model presents the relativistic 
extension of the Hartree-Fock-Bogoliubov framework, and it provides a unified description of 
ph and pp correlations. In the RHB framework the ground state of a nucleus can be written 
either in the quasiparticle basis as a product of independent quasi-particle states, or in the 
canonical basis as a highly correlated BCS-state. By definition, the canonical basis diago- 
nalizes the density matrix and it is always localized. It describes both the bound states and 
the positive-energy single-particle continuum. The matrix equations of the QRPA take a 
particularly simple form in the canonical basis, including only the matrix elements VJvk'a 
the residual p/i-interaction, and the matrix elements V^^,^y of the pairing pp-interaction, as 
well as certain combinations of the occupation factors m^, f k of the canonical single-nucleon 
states. The RQRPA configuration space includes the Dirac sea of negative energy states. 
In addition to the configurations built from two-quasiparticle states of positive energy, the 
RQRPA configuration space must also contain pair-configurations formed from the fully or 
partially occupied states of positive energy and the empty negative-energy states from the 
Dirac sea. The inclusion of configurations built from occupied positive-energy states and 
empty negative-energy states is essential for current conservation and the decoupling of spu- 
rious states. In recent applications of the relativistic (Q)RPA it has been shown that the 
fully consistent inclusion of the Dirac sea of negative energy states in the R(Q)RPA con- 
figuration space is essential for a quantitative comparison with the experimental excitation 
energies of giant resonances 0, [3, ll^ ■ 

The RHB-I- RQRPA model is fully self-consistent. For the interaction in the particle- hole 
channel effective Lagrangians with nonlinear meson self-interactions or density-dependent 
meson-nucleon couplings are used, and pairing correlations are described by the pairing part 
of the finite range Gogny interaction. Both in the ph and pp channels, the same interactions 
are used in the RHB equations that determine the canonical quasiparticle basis, and in 
the matrix equations of the RQRPA. This is very important, because the energy weighted 
sum rules are only satisfied if the pairing interaction is consistently included both in the 



5 



static RHB and in the dynamical RQRPA calculations. In both channels the same strength 
parameters of the interactions are used in the RHB and RQRPA calculations. 

In Ref. we have derived the matrix equations of the corresponding proton-neutron 
relativistic QRPA (PN-RQRPA) for an effective Lagrangian characterized by density- 
dependent meson-nucleon couplings. The PN-RQRPA has been applied in a study of charge- 
exchange modes: isobaric analog resonances and Gamow- Teller resonances. The model in- 
cludes both the T = 1 and T = pairing channels, and presents a relativistic extension 
of the fully-consistent proton-neutron QRPA that has been formulated in Ref. j^, and em- 
ployed in an analysis of /3-decay rates of r-process nuclei. In order to build a complete basis, 
the model space also includes configurations built from occu pied positive-energy states and 
empty negative-energy states. It has been shown recently jlllllJl. that the total GT strength 
in the nucleon sector is reduced by ~ 12% in nuclear matter, and by ^ 6% in finite nuclei 
when compared to the Ikeda sum rule 

{S^T-S^^)=3{N-Z), (1) 

where S^± denotes the total sum of Gamow- Teller strength for the transition. The 
reduction has been attributed to the effect of Dirac sea negative-energy states, i.e. the 
missing part of the sum rule is taken by configurations formed from occupied states in the 
Fermi sea and empty negative-energy states in the Dirac sea. 

We consider transitions between the O"*" ground state of a spherical even-even parent 
nucleus and the l"*" excited state of the corresponding odd-odd daughter nucleus, which are 
induced by the Gamow- Teller operator err . The matrix equations of the PN-RQRPA read 

(.■-;)(")-(::)(;:' 

The matrices A and B are defined in the canonical basis [16| 



A , , 

-'^prij'p'n' 



pp 



+ {UpVnUp'Vn' + VpUnVp'Un') Kf„^„p/ + {UpUnUp'Un' + VpVnVp'Vn') Kfwn' 



B. 



{UpVnVp'Un' + VpUnUpiVn') Kfi^„„' 



pn,p'n' 



- {UpUnVp'Vn' + VpVnUp'Un') Kmp'n' • 



(3) 



Here p, p', and n, n' denote proton and neutron quasiparticle canonical states, respectively, 
V^^ is the proton-neutron particle-hole residual interaction in the l"*" channel, and V^^ is 
the corresponding particle-particle interaction. The canonical basis diagonalizes the density 
matrix, and the occupation amplitudes Vp^n are the corresponding eigenvalues. However, the 
canonical basis does not diagonalize the Dirac single- nucleon mean- field Hamiltonian h£, or 
the pairing field A, and therefore the off-diagonal matrix elements H^l^i and Hp^, appear in 
Eq. (0): 

Hll, = {u^u^i - v^v^^)h^^, - {u^v^i + v^u^:)A^^> . (4) 

For each energy Ex, and in Eq. (jS)) denote the corresponding forward- and backward- 
going QRPA amplitudes, respectively. The total strength for the transition between the 
ground state of the even-even (N,Z) nucleus and a 1^ state of the odd-odd (N-1,Z-|-1) 
nucleus, induced by the Gamow- Teller operator, reads 

2 



< p\\aT^\\n > (^Xp^UpVn - Yp^VpUn 

pn 



(5) 



The spin-isospin-dependent interaction terms are generated by the tt- and p-meson ex- 
change. Although the direct one-pion contribution to the nuclear ground state vanishes at 
the mean-field level because of parity conservation, it has to be included in the calculation of 
spin-isospin excitations. The particle-hole residual interaction of the PN-RQRPA is derived 
from the Lagrangian density 

^r+p = -gpi^rp.ri; - I^^^.Yd.^fi, . (6) 

Vectors in isospin space are denoted by arrows, and boldface symbols will indicate vectors 
in ordinary three-dimensional space. 

The coupling between the /9-meson and the nucleon is assumed to be a vertex function of 
the baryon density. In Ref. p] it has been shown that the explicit density dependence of the 
meson-nucleon couplings introduces additional rearrangement terms in the residual two-body 
interaction of the RRPA, and that their contribution is essential for a quantitative description 
of excited states. However, since the rearrangement terms include the corresponding isoscalar 
ground-state densities, it is easy to see that they are absent in the charge exchange channel, 
and the residual two-body interaction reads 

V{ri,r2) = Tif2{Pj^)i{p-f^,)2gp[pv{ri)]gp[py{r2)]Dp{ri,r2) 
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rir2(SiVi)(S2V2)/^.(ri, r^) . (7) 



Dp(7r) denotes the meson propagator 



I ^-■mp(^)\ri-r2\ 

Dp(n) = , 7— , (8) 

47r ri — 7-2 



and 



cr 

For the p-meson couphng we adopt the functional form used in the DD-MEl density- 
dependent effective interaction |6|] 

9pipv) = 9pipsat)exp[-ap{x - 1)] , (10) 

where x = Pv/psat, and psat denotes the saturation vector nucleon density in symmetric 
nuclear matter. For the pseudovector pion-nucleon coupling we use the standard values 

= 138.0 MeV = 0.08 . (11) 

The derivative type of the pion-nucleon coupling necessitates the inclusion of the zero-range 
Landau-Migdal term, which accounts for the contact part of the nucleon-nucleon interaction 

Vs. = g'(^] m^i-^25in-r2) , (12) 

with the parameter g' adjusted to reproduce experimental data on the GTR excitation 
energies. 

In the pp-channel of the RHB model we have used a phenomenological pairing interaction, 
the pairing part of the Gogny force, 

VPP{1,2) = ^'^^^"'""'^'^'^^^i + BiP" - HiP^ - MiP'^P^), (13) 

j=l,2 

with the set DIS [l^ for the parameters pi, Wi, Bi, Hi and Mj {i = 1,2). This force has 
been very carefully adjusted to the pairing properties of finite nuclei all over the periodic 
table. In particular, the basic advantage of the Gogny force is the finite range, which 
automatically guarantees a proper cut-off in momentum space. For the T = proton- 
neutron pairing interaction in open shell nuclei we employ a similar interaction which was 
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used in the non-relativistic QRPA calculation [J| of /3-decay rates for spherical neutron-rich 
r-process waiting-point nuclei. It consists of a short-range repulsive Gaussian combined with 
a weaker longer-range attractive Gaussian 



Vu = -VoJ29,e-'-"'^''^tls=i,T=o , (14) 
i=i 

where fls=^T=o projects onto states with 5=1 and T = 0. We follow the prescription 
from Ref. |4J, and take the ranges /ii=1.2fm and /i2=0.7fm of the two Gaussians from the 
Gogny interaction (jl3p . The relative strengths gi = 1 and g2 = —2 are chosen so that the 
force is repulsive at small distances. The overall strength of the force represents the only 
remaining free parameter. We have already successfully applied this force in an analysis of 
the Gamow- Teller resonances in tin isotopes jll|. 



B. Calculation of /3-decay half-lives 

The rate for the decay of an even-even nucleus in the allowed Gamow- Teller approximation 
reads 



1/2 m m 



E<Ull^r-l|o+> 



pn 



where D = 6163.4 ± 3.8 s 18]. |0^) denotes the ground state of the parent nucleus, and 
1^^ is a state of the daughter nucleus. The sum runs over all final states with an excitation 
energy smaller than the Qp- value. In order to account for the universal quenching of the 
Gamow- Teller strength function, we use the effective weak axial nucleon coupling constant 
gA = 1, instead of gA = 1-26 The kinematic factor in Eq. (fT^ can be written as 



E.^El -ml{uj- E,yF{Z, A, E,) , (16) 

where u denotes the energy difference between the initial and the final state. The Fermi 
function F{Z,A, E^) corrects the phase-space factor for the nuclear charge and finite nuclear 
size effects 

C. The nucleon effective mass in the RMF models 

In nonrelativistic mean-field models the effective nucleon mass m* characterizes the en- 
ergy dependence of an effective local potential which is equivalent to the nonlocal and fre- 
quency dependent microscopic nuclear potential j2l|. m* represents a measure of the density 
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of single-nucleon states around the Fermi surface and, therefore, it has a pronounced effect 
on the calculated properties of ground and excited states. In the case of Skyrme-type in- 
teractions, for instance, calculation of ground-state properties and excitation energies of 
quadrupole giant resonances have shown that a realistic choice for the nucleon effective 
mass is in the interval m*/m = 0.8 ±0.1 [2^ . 

In the relativistic mean-field framework the expression "effective mass" has been used to 
denote different quantities. The quantity which is usually used to characterize an effective 
interaction, and which in the literature is most often called "the relativistic effective mass" , 
is also known as the "Dirac mass" [2^ 

mD = rn + S{r), (17) 

where m is the bare nucleon mass and S{r) denotes the scalar nucleon self-energy. The 
concept of the effective nucleon mass in the relativistic framework has been extensively 



analyzed in Refs. 23|,|2J]. Specifically, it has been pointed out that the Dirac mass should 
not be identified with the effective mass of the nonrelativistic mean-field models. Instead, 
the quantity which should be compared with the empirical effective mass derived from the 
nonrelativistic analyses of scattering and bound state data is given by 

m* /m = 1 — V/m , (18) 

where V denotes the time-like component of the vector self-energy. The Dirac mass, on 
the other hand, is determined by two factors: (i) the empirical spin-orbit splittings in finite 
nuclei, and (ii) the binding energy at the saturation density in nuclear matter. In the first 
order approximation, and assuming spherical symmetry, the spin-orbit part of the effective 
single-nucleon potential reads 

where M is specified as 

M = M-]^{V-S). (20) 

While the difference between the vector and scalar potentials determines the spin-orbit 
potential, their sum defines the effective single-nucleon potential and is determined by the 
nuclear matter binding energy at saturation density. The energy spacings between spin- 
orbit partner states in finite nuclei, and the nuclear matter binding and saturation, place the 
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following constraints on the values of the Dirac mass and the nucleon effective mass: 0.55m < 
rriD < 0.6m, 0.64m < m* < 0.67m, respectively. These values have been used in most 
standard relativistic mean- field effective interactions. First, we notice that in comparison 
to the nonrelativistic models, the relativistic nucleon effective mass has a rather low value, 
and this results in a smaller density of states around the Fermi surface. Second, the range 
of allowed values of the nucleon effective mass is very narrow in the standard relativistic 
mean-field phenomenology, and there is really no room for any significant enhancement of 



the single-nuc 
In Ref. b 



eon level densities at the Fermi surface, 
we have extended the standard relativistic mean-field model by including 
dynamical effects originating from the coupling of single-nucleon motion to collective surface 
vibrations. By using a simple linear ansatz for the energy dependence of the scalar and vector 
nucleon self-energies, we were able to describe simultaneously bulk nuclear properties and 
single-nucleon spectra in a self-consistent relativistic framework. 

As we have emphasized in the introduction, in order to be able to reproduce the data 
on /3-decay lifetimes, the description of single-particle energies around the Fermi surface has 
to be improved. In principle one could use the method of Ref. j3] with energy- dependent 
nucleon self-energies, but here we will try a simpler approach to increase the density of 
single-nucleon states without going beyond mean-field level. An increase of the effective 
mass necessitates a reduction of the vector self-energy (see Eq. fll8|) ). However, in order to 
retain the empirical value of the nuclear matter binding energy, the scalar self-energy should 
be reduced correspondingly. A serious problem arises because such an effective interaction 
would systematically underestimate the spin-orbit splittings in finite nuclei. A solution to 
this problem has been known for a long time, namely the tensor coupling of the tu-meson 
to the nucleon. This interaction enhances the effective spin-orbit potential in finite nuclei, 
but is not included in the most commonly used relativistic mean-field models. In Ref. j26| 
it was shown that the tensor coupling 

generates an additional term in the spin-orbit part of the effective nucleon potential, which 
now reads 



Vs.o.ir) 



1 .\d^V-S)+ 



/ ■ cr . (22) 



4M2 r ' ' 2MM r dr 

With the inclusion of the tensor omega-nucleon coupling it becomes possible to reproduce 
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the empirical spin-orbit splittings, even when using effective interactions with the Dirac 
mass as large as m^) ^ 0.7m. 

Our strategy is now to use the additional tensor-coupling term to generate a density- 
dependent relativistic effective interaction, with a value for the effective mass close to those 
used in nonrelativistic mean-field models. Starting from the DD-MEl interaction that we 
have used in the PN-RQRPA analysis of eha.ge-exehange .nodes and with the „..usion 
of the additional tensor omega- nucleon interaction (1^ . the parameters of the new interac- 
tion have been adjusted simultaneously to properties of nuclear matter and finite nuclei p. 
An additional constraint has been placed on the value of the nucleon effective mass. The 
modified effective interaction, denoted as DD-MEl*, exhibits the following values for the 
Dirac mass and the nucleon effective mass: = 0.67m, m* = 0.76m, respectively. These 
are the highest values for which a realistic description of nuclear matter and finite nuclei 
is still possible, i.e. the quality of the calculated nuclear matter equation of state and of 
ground-state properties of spherical nuclei is comparable to that of the DD-MEl interac- 
tion. The value of the Dirac mass is also in agreement with the results of Ref. {2^, where 
a detailed analysis was performed on the correlation between the isoscalar tensor coupling 
and the Dirac mass in successful mean-field models. Although the value of m* is still lower 
than those typically used in nonrelativistic mean-field models, this result presents a signifi- 
cant improvement over the standard DD-MEl density-dependent interaction (m/j = 0.58m, 
m* = 0.66m). 

We have used the new interaction to calculate the energy spacings between spin-orbit 
partner states in the doubly closed-shell nuclei ^^O, ^°Ca, '^^Ca, ^^^Sn, and ^°'^Pb. The results 
are shown in Tabled in comparison with those obtained using the DD-MEl interaction, and 
with the experimental data. Both interactions provide an excellent description of the spin- 
orbit splittings in finite nuclei. In order to illustrate the effect of the tensor-coupling term, in 
Fig. H we display the radial dependence of the spin-orbit term of the single-nucleon potential 
in the self-consistent solutions for the ground-state of ^^^Sn, calculated with the DD-MEl 
and DD-MEl* effective interactions. For DD-MEl*, in particular, also the contributions of 
the first and second term in Eq. (j^^ are plotted separately. We notice that even though the 
strength of the spin-orbit interaction that arises from the large scalar and vector self-energies 
(first term in Eq. ()22|) ) is significantly reduced, the tensor omega- nucleon coupling effectively 
compensates this reduction, and the resulting spin-orbit potential is even slightly stronger 
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than the one calculated with the DD-MEl interaction. Therefore, while both interactions 
produce very similar results for the spin-orbit splittings in finite nuclei, the inclusion of the 
isoscalar tensor-coupling term in DD-MEl* allows for an increase of the Dirac mass and 
effective mass. 

In Figs. 121 and 01 we display the neutron and proton single-particle levels in ^'^Ni and ^^^Sn, 
respectively, calculated with DD-MEl, DD-MEl*, and with two nonrelativistic Skyrme 
interactions that have been used in the calculation of /5-decay half-lives: SkO' ^ and 
SkSClT 18j. Both nonrelativistic interactions are characterized by large values of the nu- 
cleon effective mass: m* = 0.89 for SkO', and m* = 1.0 for SkSClT. Compared to the 
original DD-MEl interaction, the enhancement of the effective mass in DD-MEl* results in 
the increase of the density of states around the Fermi surface, similar to the spectra calcu- 
lated with the nonrelativistic interactions. In Table |nl we compare the calculated neutron 
and proton single-particle levels in ^^^Sn, with experimental data. The levels calculated 
with DD-MEl* are in .n„cK bette. a„t With data 0, tM.Uhose 0Ma.ed w.tK t.e 
original DD-MEl interaction. 

In Fig. IHwe plot the evolution of the vr 2^3/2 and 7^1 f 5/2 single-particle levels in heavy 
copper isotopes. The level structure has been studied in a recent experiment [28], and it has 
been pointed out that a possible reordering of these proton levels could have a pronounced 
effect on the shell structure and decay properties for nuclei near and beyond the doubly magic 
''^Ni. The results calculated with the DD-MEl* interaction are in excellent agreement with 
the experimental data. The model reproduces the observed lowering of the vrl/5/2 state in 
^^Cu, ''^Cu, and ^^Cu. For the heavier isotopes experimental data are not yet available. 
Nevertheless, our calculation predicts a further lowering of the 7rl/5/2 state, and already in 
'^^Cu isotope the state becomes the ground state. 



III. /3-DECAY HALF-LIVES 

In Ref. we have employed the proton-neutron RQRPA in a study of the high-energy 
part of the GT strength function, in particular the GT resonance which represents a coherent 
superposition of high-lying proton-particle - neutron-hole configurations associated with 
charge-exchange excitations from j = I + \ neutron orbitals into j = I — \ proton orbitals. 
In addition to the high-energy resonance, the GT strength function displays a concentration 
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of strength in the region of low-energy excitations. These transitions correspond to core- 
polarization (j = / ± i ^ j = / ± i), and back spin-flip (j = / — | ^ j = / + |) neutron- hole 
- proton-particle excitations. It is precisely this low-energy tail that contributes to the 
/3-decay process. 

Only the isovector channel of the effective interaction contributes to the matrix elements 
of the residual interaction in the calculation of charge-exchange excitations. Nevertheless, 
the isoscalar channel plays an important role because it determines the single-particle levels 
that enter the PN-RQRPA calculations. In the analysis of charge-exchange modes we have 
used the DD-MEl effective interaction, combined with the standard parameters for the pion- 
nucleon Lagrangian: m-n = 138 MeV, and /^/47r = 0.08. The parameter of the zero-range 
Landau-Migdal force g' = 0.55 has been adjusted to reproduce the excitation energy of the 
GT resonance in ^^'^Pb. We have shown that, although the inclusion of T = pairing in the 
residual interaction does affect both the low- and high-energy regions of the GT strength 
distribution, the position of the resonance is not sensitive to its strength On the other 
hand, T = pairing has a very strong influence on the low-lying tail of the GT strength 
distribution. 

Our flrst attempt to apply the PN-RQRPA to /3-decay processes, by employing the DD- 
MEl interaction, was not successful. We tried to reproduce the empirical half-life of the ^^Ni. 
This spherical nucleus undergoes high-energy fast /3-decay and, since it has doubly-closed 
spherical shells, there is no contribution from the T = pairing. The PN-RQRPA with the 
DD-MEl ...tera.fo,. pred.cts a half-life ^ 7 s, whichis an order of ,„ag„,tude longer 
than the experimental value T1/2 = 104 + 126 — 57 ms |22l|. The results of nonrelativistic 
QRPA calculations are much closer to the empirical half-life. For example, for the interaction 
DF3 the calculated half-life is T1/2 ~ 0.3 s whereas the Skyrme interaction SkO' predicts 
Ti/2 ~ 0.6 s 4]. Obviously, while the DD-MEl interaction provides an accurate description 
of the high-energy region of the GT strength function, it does a rather poor job for the low- 
energy tail of the distribution. The problem, as we have already emphasized, lies in the low 
density of single-proton states around the Fermi level. This has motivated the adjustment 
of the new effective interaction DD-MEl*, which has a considerably higher value of the 
nucleon effective mass, and consequently produces a higher density of single-nucleon states 
at the Fermi surface. 

The change of the density of single-nucleon states will also affect the position of the GT 



14 



resonance, and therefore the parameter of the zero-range Landau-Migdal force has to be 
readjusted. A higher density of states imphes a lowering of the GT strength distribution, 
and this means that value of g' = 0.55 used with DD-MEl, has to be increased in order 
to reproduce the empirical excitation energy of GT resonances. For DD-MEl* we have 
adjusted the new value g' = 0.62 to the position of the GT resonance in ^°^Pb. With this 
value of g' we also find an excellent agreement between the calculated and experimental 
excitation energies of the GTR for ^'^Ca, ^°Zr, and ^^^"^^"^Sn. 

With the new set of parameters (DD-MEl* effective interaction, = 138 MeV, 
f^/4:7[ = 0.08, g' = 0.62), we have recalculated the /3-decay half-life of ^^Ni. The value 
Ti/2 = 0.9 s presents a significant improvement over our first result obtained with DD-MEl, 
although it still overestimates the empirical half-life of ^'^Ni, and the values calculated with 
the nonrelativistic PN-QRPA. The most probable reason is that the value of the nucleon 
effective mass used in DD-MEl* is still below the values used in the nonrelativistic effective 
interactions 0,0]. A further increase of m* in our relativistic model is, however, not possible 
without downgrading the agreement with experimental data on ground state properties of 
finite nuclei. 

In recent years a number of experimental and theoretical studies have focused on the level 
structure and decay properties of neutron-rich nuclei in the vicinity of the doubly magic ''^Ni 
and ^^^Sn. In particular, beta-decay rates in these two regions of the periodic chart have 
been extensively investigated in the framework of the nonrelativistic PN-QRPA QQQ. In 
the next two sections we apply the newly formulated relativistic PN-QRPA in the calculation 
of /9-decay half-lives of nuclei in the regions Npa50 and Nf»82. 

A. Nss50 region 

In this region we have investigated the iron, nickel, and zinc isotopic chains, and the 
= 50 isotones. The structure of the low-energy part of the GT strength distribution 
crucially depends on the occupancy of the Z = 28 proton shell. It should be noted that, 
because of possible deformation effects in the Fe and Zn chains, these nuclei might not be as 
good as the Ni isotopes for a comparison with results of spherical PN-QRPA calculations. 

In Fig. El we display the calculated half-lives of the Fe isotopes (two holes in the 7^1/7/2 
orbit) in comparison with the available experimental data The results are obtained 
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with the DD-MEl* interaction, = 138 MeV, f^/An = 0.08, g' = 0.62, Gogny DIS T=l 
pairing. Two values for the strength parameter of the T = pairing interaction (fT^ have 
been used. Obviously, the absolute values of the calculated half-lives are very sensitive to 
the T = pairing strength. Without the inclusion of this pairing channel, in all considered 
cases the calculated half-lives are at least an order of magnitude longer than the experimental 
values. In the case of iron isotopes the pairing strength parameter Vq = 115 MeV has been 
adjusted to reproduce the half-life of ^*Fe. With this value we are able to reproduce the half- 
lives of the ^^Fe and ™Fe isotopes very accurately, whereas the lifetime of the ^^Fe isotope 
is somewhat overestimated. 

It is, however, probable that the inclusion of a strong T = pairing partially compensates 
the deficiencies of the single-particle spectra calculated with the DD-MEl* interaction. For 
the semi-magic nucleus ''^Fe, in Fig. El we plot the probabilities A/j of /5-transitions from an 
initial nuclear state z to a final state /, for different values of the T = pairing strength. The 
ground-state occupation probabilities for the single-particle levels relevant for this /5-decay 
process are listed in Tab. IIIII Since the vrl/7/2 orbit is not fully occupied, the transition 
with the highest probability is dominated (95% of the neutron-to-proton QRPA ampli- 
tude) by the back spin-flip configuration z/l/5/2 —>■ 7rl/7/2. Other transitions, with much 
smaller probabilities, correspond to the back spin-flip configuration z/2pi/2 — > 7r2p3/2, and 
the core-polarization configurations z/l/5/2 7rl/5/2, 1^^99/2 ^ '^^99/2, ^'^P3/2 7i"2p3/2, and 
^'^Pi/2 7r2pi/2. Since except 711/7/2, all proton single-particle levels listed in Tab. IIIII have 
very small occupation probabilities, the only sizeable contribution from the T = pairing 
to the RQRPA matrices comes from the vrl/7/2(z^l/5/2)^^ pair. Because of the attractive 
nature of the pairing interaction, the large diagonal matrix element fpV^'^p^ (p and n de- 
note vrl/7/2 and z/l/5/2 states, respectively) effectively reduces the sum of the quasiparticle 
energies: + H}^ = Ep + En- This means that the T = pairing compensates for the 
fact that even the DD-MEl* interaction still predicts a rather low density of states around 
the Fermi surface, i.e., the 7clff/2 and z/l/5/2 single-particle levels are still too close. The 
inclusion of the T = pairing will affect only configurations with the proton level at least 
partially occupied. This is clearly seen in Fig. IHl where we notice that only the transition 
built from the configuration z/l/5/2 vrl/7/2 is lowered in energy and enhanced with the 
increase of the strength of T = pairing, whereas the remaining transitions are unaltered. 

For the Ni isotopes the 711/7/2 orbit is completely occupied, i.e., in this case the tran- 
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sition 1^1/5/2 7rl/7/2 is blocked. The T = pairing could only have an effect on the 
'^^99/2(^^99/2)"^ configuration, because the i^l5'9/2 orbit is not fully occupied for isotopes 
below ^^Ni. However, it is true that when this orbit is almost empty (™Ni, ^^Ni) there 
is a large contribution from the pairing interaction, but at the same time there is only a 
small number of neutrons which can participate in the /3-decay process. On the other hand, 
when this orbit is almost full (^^Ni, ''^Ni), the contribution from the T = pairing becomes 
negligible. In contrast to the iron isotopic chain, we were not able to obtain a single value 
of the T = pairing strength parameter that would provide a consistent description of the 
entire chain of nickel isotopes. To reproduce the experimental half-lives, an extremely strong 
T = pairing would have to be used and this would cause the collapse of the PN-RQRPA 
calculation. Because of the closed Z = 28 proton shell, the T = pairing is not effective in 
Ni isotopes and, therefore, the calculated /5-decay half-lives, shown in Fig. [71 overestimate 
the experimental data. 

In this mass region we have also analyzed the half-lives of zinc isotopes. In Ref. 0] the 
measured half-lives of three zinc isotopes ''^Zn, ^^Zn, and ^°Zn, and of ^^Ge, were used to 
adjust the T = pairing strength. It was shown that the experimental lifetimes can all 
be reproduced with a single value of the T = pairing strength, Vq = 230 MeV. In the 
present calculation the measured half-lives are only reproduced by using a much stronger 
T = pairing interaction. The effect that T = pairing has on the /3-decay probability 
is illustrated in the example of the semi-magic nucleus ^''Zn. In Tab. IIIII the ground-state 
occupation probabilities for the relevant single-particle levels are included. It is important to 
note that, because of the T = 1 pairing, the occupation probability of the 1^1/7/2 state is not 
equal to one. Besides 711/7/2, only vrl/5/2 has a sizeable occupation probability among the 
proton states. The T = pairing interaction produces a large contribution to the RQRPA 
matrices for the following configurations: vrl/5/2(z^l/5/2)~^ and vrl/7/2(z^l/5/2)~^- Because 
the transition z/l/5/2 — > vrl/7/2 is essentially blocked, the effect of T = pairing is much 
weaker than in the case of Fe isotopes. If the parameter Vq is kept below ^ 250 MeV, the 
T = pairing has virtually no effect on the calculated half-lives. In the interval between 
Vo = MeV and Vq = 220 MeV the half-life decreases by just 10%. With a further increase 
of Vq, however, the calculated half-life displays a steep decrease. The measured half-life is 
reproduced for Vq ~ 330 MeV. The corresponding distributions of neutron-to-proton QRPA 
amplitudes for three values of the T = pairing strength (Vq = MeV, Vq = 255 MeV, 
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and Vq = 330 MeV) are shown in Tab. IIVI As one would expect, in the absence of T = 
pairing, the /3-decay process is characterized by two transitions. The lower one is dominated 
by the back spin-flip transition z/2pi/2 2^3/2, whereas the higher component represents a 
mixture of core-polarization transitions. Increasing the T = pairing strength to Vq = 255 
MeV, we note that: (i) transitions built from the back spin-flip configuration z/2pi/2 — ^ tt 2^3/2 
and the core-polarization configurations are no longer well separated, and (ii) an additional 
transition appears, built dominantly on the back spin-flip configuration z/l/5/2 711/7/2, 
which results in a sudden reduction of the calculated half-life. Of course, this would not be 
possible if vrl/7/2 was fully occupied. This was the case for the Ni isotopes, and consequently 
their half-lives could not be improved by increasing the strength of the T = pairing. Since 
the T = pairing has a strong effect on the z/l/5/2 7rl/7/2 configuration, a further increase 
of its strength lowers the energy of the corresponding transition. For Vq = 330 MeV we find 
only one transition, predominantly based on the z/l/5/2 7rl/7/2 configuration. 

Although Vo = 330 MeV reproduces the experimental half-lives of Zn isotopes, as shown 
in Fig. IHl this does not mean that one should use this particular value for the other isotopic 
chains, even in the same mass region. This is illustrated in Fig. IHl where we display the 
half-lives of the = 50 isotones for three values of the Vq parameter: Vq = MeV, Vq = 115 
MeV (used in the calculation of Fe isotopes), and Vq = 330 MeV (used for Zn isotopes). In 
contrast to the result of Ref. j^, the value Vq = 330 MeV which reproduces the half-lives 
of Zn nuclei, overestimates the half-life of ^^Ge. In order to reproduce the experimental 

n 

value Ti/2 = 4.55 ± 0.05 s [30], the pairing strength must be increased to Vq = 345 MeV. 
Both values, Vq = 330 MeV and Vq = 345 MeV, are rather close to the point at which the 
PN-RQRPA calculation collapses. This is illustrated in Fig. where we plot the half-lives 
of ^'^Zn, ^°Zn, ^^Zn, and '^^Ge as functions of the T = pairing strength parameter. The 
PN-RQRPA instability point is, for all four nuclei, located at Vq ^ 400 MeV. 

B. NRi82 region 

In the mass region around the doubly magic ^^^Sn, we have calculated the /3-decay half- 
lives of the cadmium, tin, and tellurium isotopic chains, as well as the A^ = 82 isotones. In 
Fig. ^2 we plot the calculated half-lives of the Cd isotopes. The results correspond to two 
calculations, with Vq = and Vq = 225 MeV for the strength parameter of the T = pairing. 
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As in the cases that we have considered in the previous section, the calculated half-lives are 
more than an order of magnitude too long when the T = pairing is not included. With 
the pairing strength parameter Vq = 225 MeV adjusted to reproduce the half-life of ^^°Cd, 
the PN-QRPA calculation reproduces the experimental half-lives of the Cd isotopic chain. 
With two holes in the 'n'lg9/2 orbit, the situation in the cadmium chain is similar to that 
of Fe isotopes (two holes in the nlfi/2 orbit). The /?-decay process in the Cd isotopes is 
dominated by the back spin-flip transition i'lg7/2 '^^99/2- Again, an increase of the T = 
pairing strength partially compensates for the fact that the difference between the vlgT/2 
and 7Tlgg/2 single-particle energies is too small, due to a relatively small effective mass. 

For the Sn isotopes the nlgg/2 orbit is completely occupied and the transition i'lgj/2 
7rl(79/2 is blocked. A similar problem we already encountered for the Ni isotopes. One 
would, therefore, expect that the calculated half-lives of the Sn isotopes will overestimate 
the experimental values by at least an order of magnitude, and furthermore that it will 
not be possible to improve the results by simply increasing the T = pairing strength. 
However, this turns out to be true only for ^^^Sn. Our model predicts that this isotope is 
stable against /5-decay, whereas the experimental half-life is T1/2 = 39.7 ± 0.5 s [3^. In the 
left panel of Fig. we display the calculated half-lives of the Sn isotopes, in comparison with 
the available experimental data js^. We notice that, in contrast to ^'^^Sn, the theoretical 
half-lives of the heavier Sn isotopes show a pronounced dependence on the T = pairing 
strength. For Vq = 225 MeV the calculated half-lives are only slightly longer than the 
experimental values. This is easily explained by the fact that in tin isotopes beyond ^^^Sn 
neutrons begin to occupy the i'lhg/2 single-particle level, and this enables the back spin-flip 
transition Z/I/19/2 —* 7rl/iii/2- The occupation probabilities of the 1^1/^9/2 single-particle level 
for the Sn and Te isotopes are included in Tab. |3 Because the neutron level z/l/i9/2 has 
low occupancy, the T = pairing produces a very strong effect on the 7rl/iii/2('^l/i9/2)~^ 
configuration, and reduces the calculated half-lives to the experimental values. 

This effect is further illustrated in the middle panel of Fig. El where we plot the calculated 
half-lives of the Te isotopes, in comparison with available experimental data. In this case, 
the choice Vq = 225 MeV results in half-lives that are even somewhat shorter than the 
experimental values. This is related to the fact that the calculated energy spacing between 
the i'lhg/2 and 'n'lhii/2 states in ^^^Sn is larger than the experimental value (see Table HHl . A 
slightly lower value of the T = pairing strength (Vq = 200 MeV) accurately reproduces the 
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half-lives of the Te isotopes. Finally, in the right panel of Fig. we display the ca^ 
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half-lives of the = 82 isotones, compared to the shell- model results of Ref. [3l|. This 
comparison also shows that Vq = 225 MeV presents a realistic choice for the strength of the 
T = pairing in this mass region. 



IV. CONCLUDING REMARKS AND OUTLOOK 



In this work the newly developed relativistic PN-QRPA has been for the first time ap- 
plied to Gamow- Teller /3-decays of nuclei relevant for the r-process. The matrix equa- 
tions of the QRPA are formulated in the canonical single- nucleon basis of the relativistic 
Hartree-Bogoliubov (RHB) model. The RHB+QRPA model employed in this work is fully 
self-consistent. For the interaction in the particle-hole channel effective Lagrangians with 
density-dependent meson-nucleon couplings are used, and pairing correlations are described 
by the pairing part of the finite range Gogny interaction. Both in the ph and pp chan- 
nels, the same interactions are used in the RHB equations that determine the canonical 
quasiparticle basis, and in the matrix equations of the RQRPA. The model also includes a 
proton-neutron particle-particle interaction (T = pairing), and here we have employed a 
sum of two Gaussian functions that was used in the non-relativistic QRPA calculation of 
/?-decay rates for spherical neutron-rich r-process waiting-point nuclei. 

Microscopic global predictions of weak interaction rates are very important, because 
most of the neutron- rich nuclei relevant for the r-process nucleosynthesis are not accessible 
in experiment. Calculated /3-decay half-lives depend on a detailed description of transition 
energies, shell structure, and on the choice of the residual interactions in the ph and pp 
channels. In the present analysis we have shown that standard relativistic mean-field effective 
interactions, adjusted to nuclear matter and ground-state properties of spherical nuclei, 
generally overestimate the empirical half-lives by more than an order of magnitude. The 
main reason is their low effective nucleon mass which, in the standard choice of interaction 
terms, is strongly related to the empirical energy spacings of spin-orbit partner states. Thus, 
in order to be able to reproduce measured /3-decay half-lives, we have to increase the effective 
nucleon mass of the relativistic mean-interaction used in the RHB calculation of nuclear 
ground state and in the ph channel of the QRPA residual interaction. In this work we have 
done it on the mean-field level, by including an additional isoscalar tensor-coupling term in 
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the Lagrangian, which allows for an increase of the Dirac mass and effective mass, while at 
the same time the new effective interaction reproduces the ground-state properties of finite 
nuclei, including the spin-orbit splittings. The new force has been adjusted starting from 
our most successful parameter set so far, the density-dependent interaction DD-MEl. We 
note, however, that the inclusion of the tensor term allows only for a moderate increase of 
the effective mass. We cannot, like in Skyrme forces for example, adjust interactions with 
values of the effective mass close to one. 

With the new density-dependent interaction DD-MEl* we have calculated the GT dis- 
tribution strengths and /3-decay rates of neutron-rich nuclei in the mass regions N;^50 and 
N^82. The model reproduces in detail the data on GT resonances and the low-energy GT 
strength (see also Ref. jlll]). The results for /5-decay half-lives are similar to those obtained 
in the nonrelativistic PN-QRPA calculation of Ref. ^, where it was shown that a fine tun- 
ing of the strength of the T = pairing interaction is necessary in order to reproduce the 
experimental data. In general, by adjusting the strength parameter of the proton-neutron 
pairing interaction to one experimental half-life, the PN-RQRPA calculation reproduces the 
data for a chain of isotopes. In the region N?«50 very different values of the T = pairing 
strength reproduce the empirical lifetimes of the Fe and Zn isotopic chains, whereas a single 
value, adjusted to the half-life of ^^°Cd, qualitatively reproduces the data in the N;^82 re- 
gion. T = pairing, however, does not help in the case of Ni isotopes and for ^'^^Sn, and the 
model overestimates the half-lives on Ni nuclei and predicts a /3-stable ^^^Sn. Therefore, we 
have not been able, as suggested in Ref. 0, to adjust the strength of the T = pairing on 
experimental data, and extend the calculation to other mass regions, or even to other iso- 
topic chains in the same mass region. This problem could be due, however, to the deficient 
description of shell structure, related to the low effective mass of our ph interaction. This is 
also the reason why we did not attempt a calculation of /3-decay rates in the region N;^126, 
where there are not enough data to constrain the strength of the T = pairing interaction. 

The relativistic PN-QRPA nevertheless presents a valuable tool for the investigation of 
weak interaction rates of neutron-rich nuclei. It can be used for microscopic calculations 
not only of Gamow- Teller, but also of first-forbidden decays, which play an important role 
for /?-decay rates both in the N?«82 and N^126 regions 0, In order to be able to make 
reliable predictions, however, the model will have to be improved by allowing for higher 
effective nucleon masses, possibly by going beyond the simple mean-field approximation. 
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FIG. 1: Radial dependence of the spin-orbit term of the single- nucleon potential in self-consistent 
solutions for the ground-state of ^'^^Sn. The curves correspond to RMF calculations with the DD- 
MEl and DD-MEl* effective interactions. For the latter the contributions of the first and second 
term in Eq. (|22|) are plotted separately. 



FIG. 2: Neutron and proton single-particle levels in ^^Ni calculated with DD-MEl (a), DD-MEl* 



(b), and with two nonrelativistic Skyrme interactions: SkO' ^| (c) and SkSC17 



1^ (d). 



FIG. 4: vr2p3/2 and vrl/5/2 single-particle levels in heavy copper isotopes, calculated with the 
DD-MEl* interaction, in comparison with experimental data 28]. 



FIG. 5: Comparison of the calculated half-lives of Fe isotopes, for two values of the T = pairing 



strength, with available experimental data 



FIG. 6: Probabilities of /^-transitions for the semimagic '''^Fe, for different values of the T = 
pairing strength. 



FIG. 7: Calculated half-lives of Ni isotopes, in comparison with experimental values. The data are 
from Ref. 30 1, except for ''^Ni where the value = 104 -|- 126 — 57 ms [2^ is used. 



FIG. 8: Calculated half-lives of Zn isotopes, without anc 



comparison with experimental and extrapolated values 



3C|. 



with the inclusion of T = pairing, in 



FIG. 9: Calculated half-lives of = 50 isotones for three values of the T = pairing strength. 

30|, except for '^^Ni where the value = 104 -|- 126 — 57 ms 29] is used. 



The data are from Ref. 



FIG. 10: Calculated half-lives of ^^Zn, ^°Zn, ^^Zn, and ^^Ge as functions of the T = pairing 
strength. 



FIG. 11: Calculated half-lives of Cd isotopes for two values values of the T = pairing strength, 
in comparison with available experimental data jso! ]. 



FIG. 3: Same as in Fig. HI for ^^"^Sn. 
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FIG. 12: Calculated half-lives of Sn (left panel) and Te (middle panel) isotopes for two values of 
the T = pairing strength, in comparison with experimental data y. In the right panel the 
results for the N = 82 isotones are compared with the shell- model results [si]. 



TABLE I: Energy separation (in MeV) between spin-orbit partner states in doubly closed-shell 
nuclei, calculated with the DD-MEl and DD-MEl* interactions, and compared with experimental 



data 
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DD-MEl 


DD-MEl* 


Exp. 


160 


ulp 


6.32 


6.02 


6.18 




Trip 


6.25 


5.96 


6.32 


40Ca 


vld 


6.57 


6.59 


6.00 




Tild 


6.51 


6.51 


6.00 


48Ca 




7.69 


7.79 


8.38 




v2d 


1.72 


0.56 


2.02 


i32Sn 


v2d 


1.88 


2.03 


1.65 




■Klg 


6.24 


6.57 


6.08 




TT2d 


1.82 


1.98 


1.75 


208pb 


u2f 


2.20 


2.38 


1.77 




uli 


6.84 


7.13 


5.84 






0.88 


0.89 


0.90 




TT2d 


1.65 


1.82 


1.33 




vrl/i 


5.84 


6.06 


5.56 
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TABLE II: Neutron (left panel) and proton (right panel) single-particle energies Enij in I32sn. 
For each panel, the left-hand column specifies the radial, orbital and total angular momentum 
quantum numbers, the column labeled DD-MEl contains results calculated with the standard 
density-dependent RMF interaction DD-MEl, the spectrum calculated with the new effective in- 
teraction DD-MEl* is included in the third column. The theoretical spectra are compared to the 
experimental energies shown in the column labeled EXP. 





NEUTRON STATES 






PROTON STATES 




nlj 


DD-MEl 


DD-MEl* 


EXP 


nlj 


DD-MEl 


DD-MEl* 


EXP 


I59/2 


-18.76 


-17.44 




2p3/2 


-18.91 


-18.79 




2c?5/2 


-10.91 


-10.32 


-9.04 


2pi/2 


-17.62 


-17.44 


-16.13 


157/2 


-12.86 


-11.20 


-9.82 


1^9/2 


-16.25 


-16.55 


-15.78 


3si/2 


-8.68 


-8.12 


-7.72 


I57/2 


-10.04 


-9.98 


-9.65 


1^11/2 


-7.42 


-7.30 


-7.63 


2(^5/2 


-7.60 


-8.58 


-8.69 


2(^3/2 


-8.96 


-8.29 


-7.39 


2c^3/2 


-5.70 


-6.60 


-6.95 


2/7/2 


-1.40 


-1.29 


-2.45 


1^11/2 


-5.30 


-6.78 


-6.86 


3p3/2 


-0.58 


-0.35 


-1.59 


3Sl/2 


-4.88 


-5.87 




1^9/2 


-0.11 


0.45 


-0.88 


2/7/2 


2.66 


1.32 




3pi/2 


-0.16 


0.05 


-0.75 










2/5/2 


0.38 


0.56 


-0.44 
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TABLE III: Occupation probabilities of neutron and proton single-particle states for the ground- 
states of '^^Fe and ^°Zn. 





76Fg 






80Zn 




nlj 


NEUTRONS 


PROTONS 


nlj 


NEUTRONS 


PROTONS 


1/7/2 


1.000 


0.743 


1/7/2 


1.000 


0.980 


1/5/2 


1.000 


0.019 


1/5/2 


1.000 


0.260 


2p3/2 


1.000 


0.005 


2p3/2 


1.000 


0.098 


2pi/2 


1.000 


0.004 


2Pl/2 


1.000 


0.040 


1^9/2 


1.000 


0.004 


159/2 


1.000 


0.014 



TABLE IV: /? -transitions in Zn for three values of the T = pairing strength. The contribution 
(in percent) of a particular configuration to the QRPA amplitude is included in parenthesis. 



Vo = MeV (Ti/2 = 


-- 18.9 s) 


Vo = 255 MeV (T1/2 = 13.1 s) 


Vo = 330 MeV (T1/2 


= 0.6 s) 


Ef-Ei (MeV) 




Ef - Ei (MeV) 


Ef - Ei (MeV) 




-3.088 gg/2 ^ 


59/2 (12%) 


-3.271 /5/2 ^ /5/2 (11%) 


-5.273 /7/2 ^ 


• /7/2 (5%) 


Pl/2 


P3/2 (80%) 


59/2 59/2 (18%) 


/5/2 ~^ 


/7/2 (71%) 






Pl/2 P3/2 (60%) 






-2.690 P3/2 - 


. P3/2 (6%) 


-2.982 /5/2 ^ h,2 (13%) 






Pl/2 


Pm (11%) 


59/2 59/2 (16%) 






fb/2 


/5/2 (22%) 


Pl/2 ^ P3/2 (19%) 






59/2 


59/2 (54%) 


P3/2 ^ /5/2 (47%) 










-2.302 Pl/2 ^ Pl/2 (13%) 










59/2 59/2 (15%) 










/5/2 ^ /7/2 (58%) 
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TABLE V: Occupation probabilities of the i'lhg/2 single-particle state for the ground-states of the 
Sn and Te isotopes. 





Sn 




Te 


A 




A 


ulhg/2 


134 


0.024 


136 


0.035 


136 


0.054 


138 


0.077 


138 


0.092 


140 


0.129 


142 


0.142 


142 


0.192 


144 


0.203 


144 


0.263 


146 


0.271 


146 


0.339 
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